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We investigate the behaviour of a model with two oppositely charged scalar fields. In the Bo- 
gomol’nyi limit this may be seen as the scalar sector of A^=l supersymmetric QED, and it has 
been shown that cosmic strings form. We examine numerically the model out of the Bogomol’nyi 
limit, and show that this remains the case. We then add supersymmetry-breaking mass terms to 
the supersymmetric model, and show that strings still survive. Finally we consider the extension 
to N=2 supersymmetry with supersymmetry-breaking mass terms, and show that this leads to the 
formation of stable cosmic strings, unlike in the unbroken case. 


I. INTRODUCTION 

Supersymmetric (SUSY) theories often have a continuous degeneracy of inequivalent ground states. Classically this 
corresponds to flat directions of the scalar potential, which may arise from abelian theories when the gauge symmetry 
is broken by a Fayet-Iliopoulos D-term. Recently, these have undergone intense study, mostly in the context of effective 
actions for SUSY non-abelian gauge theories and confinement [!]. 

We naturally expect the formation of cosmic strings when the D-term breaks the gauge symmetry [2, 3], but the 
presence of the flat directions means that the situation is more complicated. It is not immediately obvious which 
vacuum states among the degenerate vacua will give rise to stable string solutions. 

The situation has been investigated in [2] where a model with two oppositely charged chiral multiplets and zero 
superpotential is discussed. They show that one particular vacuum is favoured, and that a Nielsen-Olesen (NO) [4] 
type cosmic string can form in one of the scalar fields, the other field being zero everywhere. It was later shown in 
[5] that the crucial role in this vacuum selection effect is that played by the gauge field. The vaccum is effectively 
chosen so that the vector mass is minimized, leading to a lower energy for the core of the vortex. The authors in [5] 
also proved that the preferred cosmic string was stable, using a Bogomol’nyi type of argument. 

In this paper we consider the scalar sector of this theory outside the Bogomol’nyi bound to study whether the 
result mentioned above holds as well for other values of the parameters. Although it no longer corresponds to the 
simplest SUSY model, it still possesses flat directions, and may be viewed as illustrating the behaviour of more general 
SUSY theories with flat directions. An analytical approach is not possible out of the Bogomol’nyi bound, so we carry 
out numerical simulations of the fields to see if there is still any evidence of vacuum selection occurring. In fact the 
simulations suggest that, for all the range of parameters considered, the system chooses one out of all the possible 
minima of the flat directions and forms NO strings. 

We then study the fate of the different defects formed in the A^=l and N=2 supersymmetric QED models considered 
in [5] in which we break supersymmetry by the addition of soft supersymmetry-breaking mass terms to the scalars. 
Then, the degeneracy in the minima of the potential is lifted, and we have no flat directions. For the iV=l case, the 
NO string which is present without mass terms can survive SUSY breaking for some values of the masses. 

The situation is different for the N=2 model. The system was first considered in [6], in which it was claimed that 
the model accommodated strings that could confine magnetic charges. However, [7] showed that the strings were 
not stable, and so could not confine. A further step attemping to get stable strings by adding a Fayet-Iliopoulos 
term [5] proved to be insufficient. Although there is a vacuum selection effect in the model, it is not enough to 
form stable cosmic strings. The strings formed are semilocal [8], and are only neutrally stable, in the sense that 
small perturbations tend to cause the string core to expand indefinitely [9]. We will show that with the additional 
supersymmetry breaking mass terms, stable cosmic strings are possible. 


II. A SCALAR MODEL OUT OF THE BOGOMOL’NYI BOUND 

The model we consider has two complex scalar fields </>+ and </>_ of opposite charges, coupled to a U{1) vector field. 
The Lagrangian is of the form 


C = ^ , 


( 1 ) 
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where is the U{1) gauge held and, as usual, 

D^(j)+ = {d^ + iqA^) (j>+ ; 

= {df, - iqAf,) (j)_ ; 

= d^A^ - d^A^ . ( 2 ) 

The above can be regarded, in the Bogomornyi limit A = q^, as the bosonic part of a four-dimensional model 
with two fV=l chiral superhelds (d)_|., $_) of opposite charges ±g coupled to a U{1) vector multiplet, with zero 
superpotential and a Fayet-Iliopoulos term proportional to if'. 

The vacuum manifold of this model is given by 

= (3) 

with solutions 

|(/)+| = Tycoshrt = v+, |^_| = ? 7 sinhu = w_ , (4) 

where u parametrises the moduli space. When u = 0, the held is zero, and formation of the usual NO string [4] 
in the held may occur, since we are back at the abelian Higgs model. 

By considering the energy per unit length of straight static vortices, we hnd that the conditions for hnite energy 
impose some relations between the helds: since we are looking to construct cylindrically symmetric conhgurations, 
helds with charge ±q wind as and the gauge held tends to a constant Ag —n. Bearing this in mind, we use 

the ansatz 

(/.+=/+(r)e*«"^ (/._ =/_(r)e-*«"''e*^, Agd0 = a{r)Ae, (5) 

where f± and a are real functions, and A is a real constant. /±(0) = a(0) = 0 for the functions to be well-dehned, 

and a(r ^ oo) = —n by the previous energy arguments. Such considerations also tell us that the scalar helds should 
go to a minimum of the potential energy as r ^ oo, and hence the boundary conditions for the f± functions are 
f±{r oo) = u±. It is possible to construct vortices tending to any of the vacua (4) as r ^ oo. 

However, it has been shown [2, 5] that when the Bogomol’nyi limit of this system (/3 = \/q^ = 1) is considered, 

the only static solutions are those with u = 0. After symmetry breaking, apart from the massless helds, we are left 
with a massive vector held (m^ = 2 g^? 7 ^ cosh 2u) and a massive scalar held (m^ = 2Xr]^ cosh 2m). The system chooses 
M = 0 in order to minimize the vector mass. 

Thus NO strings may form only in the (p+ held, since f-{r) = 0. In the following section we would like to move 
out of the Bogomol’nyi limit to see whether m = 0 is still preferred by the system. In fact, as explained in [2], the 
dynamics of the magnetic core and those of the scalar helds outside the core are effectively decoupled. Far from the 
core, the magnetic held is very small, and the scalar helds lie on the vacuum manifold. Hence, scalar helds are free to 
move along the moduli space with no appreciable cost in energy, so one might presume that some (jj- lumps would 
be formed far from the core of the string and this might have some cosmological consequences. 


III. NUMERICAL SIMULATIONS 

We use numerical simulations of a discretized version of the model described by (1) in order to study its behaviour 
out of the Bogomol’nyi bound. We hrstly set 77 = g = 1 by rescaling the helds and spacetime coordinates, and then 
we use techniques borrowed from Hamiltonian lattice gauge theory (see for instance [10]) to perform the simulations. 
The lattice link and plaquette operators are dehned as 

( 6 ) 

Qv = Uj{x)Ui{x + Xj)Uj{x + Xi)U}{x) , (7) 

respectively, where I is the lattice spacing, the label i takes the values 1,2,3 corresponding to the three spatial 
dimensions, and Ai are the gauge helds. By x + Xi, we denote the nearest lattice point in the i direction from x. The 
plaquette operators are related to the gauge held strength. The lattice link operator, on the other hand, is used to 
dehne discrete covariant derivatives 

Di(j)+{x) = j{ul{x)(j)+{xFxi) - (l)+{x)^ ; 

Di(j)-{x) = j{Ui{x)(j)-{x + Xi) - (j)-{x)) . 


( 8 ) 
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We can then obtain a discretized Hamiltonian density corresponding to the Lagrangian (1) 




+|A</>+|^ 


\D,6-\ 




2/4 


i¥=j 


( 9 ) 


where n± and i?* are the conjugate momenta of (/)± and Ai respectively, and l3=X/q^ (/3=A in rescaled units). In the 
I —> 0 limit, the continuum Hamiltonian for the system is recovered. 

One interesting property of (9) is that it is gauge invariant, in the sense that a general C/(l) transformation A(a;) 
for which the fields transform as 


(j)+{x) - 

A(x)V+(a;); 

(j)-{x) - 

K{x)(j)-{x ); 

Ui{x) - 

K(x)Ui{x)K^{x + Xi) 


( 10 ) 


leaves the Hamiltonian invariant. 

In order to get the Hamiltonian in the form (9) a gauge choice has been made, namely, Aq = 0. The equation of 
motion corresponding to Aq is equivalent to Gauss’s law 


^<)'+(a;)n_|_ — n)^(/)_i_(a;) + n^(a;)(/()_(a;) — (}()^)(x)n_(x)^ = (x) — {x — Xi)) . 

i 


( 11 ) 


The system is geometrically forced to preserve Gauss’s law, and, if the initial conditions satisfy relation (11), subse¬ 
quent evolution by the equations of motion coming from (9) will also satisfy (11). 

The gauge invariant dissipative terms yHj- and jE^ were added to the Hamiltonian equations of motion for (j3± and 
Ai respectively. These obey Gauss’s law, and thus the system as a whole will still obey Gauss’s law as well. 7 is a 
free parameter in the model, but simulations using different values of 7 show qualitatively similar behaviour. 

The equations of motion coming from the Hamiltonian (9) were simulated numerically on a 64^ lattice. We are 
interested in the role played by the field in the dynamics of the system, not in the details of the phase transition 
itself. Furthermore, previous work concerning the numerical simulation of defects showed that the formation of 
defects depends on the later interplay between scalar and gauge fields, and is not actually very sensitive to the initial 
configuration used [11, 12]. We did, however, use different initial conditions to verify that the result was not dependent 
on how the system was set. In fact, all simulations using different initial conditions give rise to qualitatively the same 
results. The first time steps are a transient period during which the system loses energy rapidly, and relaxes. Then, 
the interaction between gauge and scalar fields leads to the formation of defects. 

All the results shown in this paper are obtained by using (f>± = 0, = 0 and A^ = 0 and random initial velocities 

for the scalar helds as the initial configuration. This choice of conditions obeys Gauss’s law, and hence we have 
guaranteed that Gauss’s law holds throughout the simulation. Gomputation of Gauss’s law during the evolution of 
different initial conditions was used to ensure the stability of the code. 


IV. NUMERICAL RESULTS 

Figure 1 is a schematic example of the outcome of the simulations. It can be seen that there exists a transitional 
period during the early timesteps, when the system is trying to get a physical configuration by dissipating energy. 
Both scalar fields acquire non-zero values at this time, and in particular Some of the initial kinetic energy 

is transformed to potential energy, but there is no hint of structure of any kind yet. At the end of the transitional 
period, \(j)-\ decreases rapidly until 1 (/()_|~ 0 . 

In a) we have plotted the number of sites in the lattice with |(/)+|<0.75, and also the number with |(()_|>0.1, as a 
fraction of the total number of lattice sites. After the transitional period (f~15), there are no points on the lattice 
with |(/)_|>0.1, whereas the field (j)+ has modulus greater than 0.75 almost everywhere. The scalars fields should tend 
to minimize the potential energy, and according to these results, it seems that only (/)_|_ is involved in achieving this. 
This is strongly suggestive of the idea that, out of all the possibilities, the system has chosen u = 0 as the preferred 
point in the flat directions (4). 

In figure I b), we have plotted the quantity 1 — where V{x) is the potential energy at lattice point x, 

and Vjnax is defined as the total potential energy calculated in the case when (j)± = 0 everywhere. |(/)+(x)| and 
\4’-{^)\ are also plotted as fractions of the total number of lattice points. In this hgure we see that, as in a). 
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after the transitional period (j)- tends quickly towards zero, while \4>+ \ ~ 1 almost everywhere. The potential energy 
mimics the behaviour of 0_|_, and ignores the presence of (/)_, suggesting once again that the actual minimization of 
the potential energy is due to the field (/>+ alone. 




FIG. 1: Representation of the results of a numerical simulation performed in a 64® lattice, with /3=0.3 and 7 = 0 . 5 . Figure a) 
shows the number of lattice sites with j</)+|<0.75 (dashed line) and |<()_1>0.1 (continuous line) as a fraction of the total number 
of sites. Figure b) shows the sum of the modulus of each of the fields at every site (as a fraction of the total lattice points); 
50^ \(f>+ \ (dashed line), \<j)- \ (continuous line). There is also a representation of the quantity 1 — V(x) /Rmax, where V(x) 

is the potential energy (dotted line). 


The system is not exactly at the minimum of the potential everywhere, as there are a number of points in the lattice 
at which |(()+|<0.75. This indicates that there are points at the false vacuum ((()+=(/)_=0), and there may be strings 

forming in the system. In figure 2 we give a representation of the magnetic field, ^, and the modulus of the 
field (/)_!_ at time t=40, after the transient period is over. It is clear that the points where |(()+|<0.75 form structures, 
namely strings, and that the areas in which the magnetic field is concentrated corresponds exactly to these points. 
A similar plot using |0_| would show that |0_|~O everywhere. The excellent match described above between the 
magnetic field and \4>+\ is another good check of the stability of our code. We performed several simulations varying 
(3, from (3=0.1 to (3=2.0, which showed the same behaviour. 



FIG. 2: Isosurface representation of the modulus of (j>+ (left), where \4>+\ ^ 0.75 is shown, and 
^ 0.25 of the measured maximum is displayed, for a 64® simulation with /? = 0.3, 7 = 


magnetic field (right), where 
0.5 at t = 40. 


These results show that after some initial (non-physical) period of relaxation, only one of the scalars becomes 
dynamical, and that the other (c3_) goes rapidly to zero everywhere, decoupling from the other fields. Thus the 
analytical result obtained for (3=1 [5] has been shown to be valid for cases with In all cases the system chooses 
the vacuum with m= 0, and the field (j)- is effectively decoupled from the dynamics. After this decoupling, we are 
left with what essentially is an abelian Higgs model and therefore the system accommodates stable Nielsen-Olesen 
strings, formed by interaction of the gauge field with 
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V. ADDING MASS TERMS 

In the previous section we demonstrated that although the vacuum of the system investigated has fiat directions, 
the only vacuum configuration appearing in the evolution of a string network is the one with m=0 , which is a stable 
cosmic string. We will now proceed to show another situation in which stable NO strings form in this context. We 
first add mass terms to our Lagrangian (1), obtaining NO vortices again, and then move on to a more complicated 
system with N=2 SUSY, which we will again break with mass terms. In this second case, when there are no mass 
terms, the strings formed are semilocal strings [ 8 ], which in the limit considered are neutrally stable. There is a family 
of such strings, all degenerate in energy, with different widths ranging from that of the NO vortex to an infinite core 
[14]. There exists a zero mode corresponding to changes in the width of the core, which when excited causes any 
given string in this family to grow indefinitely. 

Since SUSY is clearly not manifest at low energies, it is necessary that it be broken at some scale. Adding these 
soft supersymmetry-breaking mass terms to the Lagrangian is a convenient way of doing so while avoiding the details 
of how the breaking actually occurred. 

The new Lagrangian, modified from (1), is 

£ = \D^(l)+\^ + ^ - TO+|(/)+P - , (12) 

and the two extrema of the potential are now at 

((.+ =0, (13) 


and 


<(._ = 0. (14) 

When \rj^>m\, the second case is clearly the minimum. The degeneracy of the vacuum is lifted, m = 0 is selected 
out of the possible vacua (4), and we are left with a genuine NO vortex. If the condition does not hold, the 

only solution has both fields equal to zero, and no strings form. Therefore, depending in the value of the system 
either accomodates a NO string solution in (/)_|_, or no string solution is possible at all. As in the previous case, (j)- 
does not play any role in the defect formation. 

A more interesting result is obtained when this A^=l model is upgraded to N=2. We consider a model consisting of 
two N=2 hypermultiplets ha, where a=l,2, of opposite charges qa=Fq, coupled to an N=2 abelian vector multiplet. 
Such a model was considered in [ 6 ] as a candidate for obtaining magnetic flux confinement in the low energy limit 
of type II superstrings compactified on Calabi-Yau manifolds. The vortex solution obtained in [ 6 ] was found to be 
unstable [13], since the U(l) gauge symmetry was not broken. A new attempt to obtain stable vortices in this model 
[5] added a Fayet-Iliopoulos D-term of the form k-D to the system. This gives rise to a term that spontaneously breaks 
the gauge symmetry. The defects obtained in this model were the so-called semilocal strings [ 8 ], a one-parameter set 
of magnetic ’’vortices”, which includes the NO vortex. In the limit studied, semilocal strings are known to spread out 
to greater and greater radius, as mentioned before, until the string relaxes to the vacuum [14], although the magnetic 
field remains quantized. 

In the present work, we add soft supersymmetry-breaking mass terms to the Higgs scalars of the model studied in 

[5], and investigate the resulting vortices. 

-♦ 2 

Without loss of generality we take fc = (0,0, ^). The energy of the system is 


E = Jd^x^[\D^hn\^ + \D^hi2\^ + \D^h2i\^ + \D^h22\^ + 
+ - Hi + ^)2 + {Hi + Hif + (* Hi - i Hif 

Fmli\hii\^ + mi2|/il2p + TO2 iI^2iP + TO22|/i22p , 


(15) 


where D^ = d^ + i{qa)A^ and * = -{qa/‘F) hl^haj- 

The original model had, in addition to the U{1) gauge symmetry, an effective global SU{2) symmetry of the 
scalar potential, mixing the pairs of fields {hn and / 122 ) (^12 and / 122 ) amongst themselves. This symmetry was 

responsible for the appearance of semilocal strings. It would therefore be expected that the addition of different mass 
terms mu 7 ^ m 22 would lead to the formation of NO vortices, as the SU{2) symmetry is broken by those terms. 



6 


Parametrising the scalars as hai = raie^^°'', the potential takes the form 

^ ^ + rh + r?2 - ^22 + + 4rnr^2 

8 (riiri2r2ir22COs(xii + X22 - X12 - X21) ) ) 

|l/' 22 i 22 , 22 , 22 \ 

+ 2 + "il2^’l2 + W2ir2i + TO22r22) • 


+ 4riir^2- 


(16) 


In order to obtain the minimum potential we require cos(xii + X 22 — X 12 ~ X 21 ) is maximised, so that xii + X 22 — 
X 12 ~ X 21 = 0 modulo 27r. The vacuum manifold can then be calculated by varying V over all the fields, = 0, 
and we find the following minima 

2 2 o 2 

rii = 77 -2mii, 

^12 ~ ~ 2 m .22 , 

2 2 o 2 

^21 ~ V 2 m 2 i, 

2 2 o 2 

^22 ~ V 2 m 22 , 

together with the maximum 

rii = ri2 = r2i = r22 = 0, (18) 

supposing that all masses are different. We see that solutions with ri 2 7 ^ 0 and r 2 i 7 ^ 0 are never valid, so we are left 
with only two possible minima. These two minima correspond to NO vortices in rn or r 22 , and depending on the 
values of mu and TO 12 one or the other will form. 

For mil < TO 22 and 2mh < rf' the NO strings will form in the hn field; on the other hand, if TO 22 < mu and 
2 to22 < 'rf' the NO will be in the /122 field. If the masses are bigger than rf, then no vortex solution is possible. 

A special case is when mu = TO 22 < rf (2, since in this instance the SU(2) symmetry between hn and ^,22 is not 
broken, and we recover semilocal strings as in [5]. The minima in this case are 

+ ^22 =v‘^ - 2TOn , (19) 

and we have a whole family of solutions [14], the NO string among them. 

We then conclude that, depending on the values of the added masses mu and m 22 , we might have the case where 
no strings at all form, the case where we recover semilocal strings, or the system might be able to form strings in hn 
or /122 that are stable NO vortices. 


?"i2 = r2i = r22 = 0 ; 

Til = r2i = r22 = 0; 

Til = ri2 = r 22 = 0 ; 

rii = ri 2 = r 2 i = 0, (17) 


VI. DISCUSSION 

We have seen that, for a general class of scalar fields with fiat directions in the potential, there does seem to be a 
vacuum selection effect taking place even when the system is out of the Bogomol’nyi limit. Furthermore, the vacuum 
chosen remains the same as in the Bogomol’nyi limit, regardless of the value of the parameter f. After spontaneous 
symmetry breaking, defects may form. Due to this vacuum selection effect only one of the fields becomes dynamical 
and this field gives rise to Nielsen-Olesen strings, whereas the other field decouples. This is what is seen in simulations. 

On the other hand, for supersymmetric models with flat directions, the addition of supersymmetry-breaking mass 
terms, that arise in the effective SUSY action at low energies, lifts these flat directions. In the case that mii=m 22 , 
the string is a semilocal string in the Bogomol’nyi limit. This solution has the same form as that of the model without 
masses, where there are no stable string solutions. Nevertheless, if the masses are different, only one scalar field has 
a non-zero vacuum expectation value. This field then gives rise to a NO cosmic string. 

In addition, in supersymmetric theories there will be fermion zero modes in the string cores due to the symmetry 
between bosons and fermions [3]. In the case where the supersymmetry is broken, it may still be possible to have 
fermionic zero modes on these strings. The zero modes give rise to currents that can stabilise string loops. We will 
investigate this possibility in a subsequent paper [15] . 
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